arXiv: 1509.05418vl [hep-th] 17 Sep 2015 


Prepared for submission to JHEP 


A breakdown of a universal hydrodynamic relation in 
Gauss-Bonnet gravity 


Evgeny Shaverin 

Department of Physics, Technion, Haifa 32000, Israel 

E-mail: evgenyStx. technion .ac.il 

Abstract: We compute the second order transport coefficients of a hydrodynamic theory 
with Einstein-Gauss-Bonnet gravity dual. We show the breakdown of the universal hydrody¬ 
namic relation —2Ao + 4Ai — A 2 = 0 for the general values of the Gauss-Bonnet coupling. 
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1 Introduction 

Conformal viscous hydrodynamics 

Hydrodynamics is a low energy effective description of many-body systems which are studied 
at length scales much larger than their mean free path. When only energy and momentum 
are conserved, this effective description is given by a handful of hydrodynamics fields: the 
velocity field (normalized so that = — 1), and the temperature field T. The dynamics 
of these fields is governed by energy conservation, 

df,Tf^’^ = 0. ( 1 . 1 ) 

In a conformal theory, the dependence of the stress-energy tensor on the hydrody¬ 
namic fields takes the following form 

= P (4u^m^ + <^) + , (1.2) 

where P = P{T) is the pressure and H^^ denotes dissipative terms. The tensor H^'^ depends 
on gradients of the hydrodynamic variables T and u^. It is convenient to impose on H^^ the 
Landau frame condition = 0. Following [1], an expansion of H^^ in derivatives of T 

and up to second order is given by 

4 

i=0 
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(1.3) 




The coefficient rj in (1.3) is the shear viscosity, and the coefficients Aj, z = 0,..., 4, are referred 
to as ‘second order transport coefficients’. The tensors and are given by: 

S® = , 

u) 5 v) •> 

with the following definitions: R^y and Ra^vp are the Ricci and Riemann tensor respectively, 
the vorticity, Uf^y, is given by = \P^°'Py^ {daUp — dpUa), and triangular brackets denote 
a traceless and transverse projection 

Aif^y) = ^p^apup ^ Ap^) - , 

where we have defined the projection on the spatial directions + u^u’'. We assume 

that the field theory is living in flat space, where the Ricci and Riemann tensors vanish, 
therefore = 0. Collecting (1.2) and (1.3) we find that in our setup. 


3 

T^^‘' = P (4n^u^ + + Y1 ■ 

i=0 


(1.5) 


In [2] it was conjectured that for all known physical systems the ratio of the shear viscosity 
to entropy density is given by 

( 1 - 6 ) 

s dvr 

The bound (1.6) is found to be saturated in all thermally equilibrated systems with isotropic 
and homogeneous configuration described by a dual two-derivative gravity action [3-5]. The 
bound (1.6) has been shown to be violated in the non-isotropic configurations, see [6-11]. 
Recent investigations [12-14] show that in theories dual to higher derivative gravity there is 
a violation of this bound. In the special case of Gauss-Bonnet gravity, 


^egb = — 


1 


IbvrGs 




_, ^ 12 

—g i? -|- --T7 — 


L2 


eicBj 


(1.7) 


LgB = Rmnp,R^^^^ - ^RmnR^^ + R^ , 

it was shown in [13] that for general values of a Gauss-Bonnet coupling 9 

2 = 2 -( 1 - 89 ). ( 1 . 8 ) 

s dvr 

When considering the second order transport coefficients, it has been observed that the 
linear relation 

-2Ao + 4Ai-A2 = 0, (1.9) 
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seems to be satisfied for a large class of two derivative gravity theories [15]. In contrast, the 
coefficients Aj where also computed perturbatively using kinetic theory where the relation 
(1.9) seems not to hold. 

Surprisingly, Gauss-Bonnet gravity seems not to correct the relation (1.9) when it is 
treated perturbatively [16]. This result has been conhrmed and extended in [17] to include 
perturbatively small eight derivative corrections. In what follows we extend these results and 
compute the expression (1.9) to all orders in the Gauss-Bonnet coupling. We find that 


— 2Ao -|- 4Ai — A 2 


TT^r^ V2 {1 - y/1 - 89 + 169) 6 

[ (l + Vl- 80 )^^^ 


( 1 . 10 ) 


Thus, the relation (1.9) which seems to hold even at the perturbative level in 2, 4 and 8 
derivative gravity breaks down at the non perturbative level. We note that the relation (1.10) 
has been published before in [17] using methods and the non-perturbative expressions for all 
Gauss-Bonnet transport coefficients published in [18]. Treating the Gauss-Bonnet coupling 
9 non perturbatively via AdS/CFT may be regarded as a toy model for probing finite N 
corrections in a gauge theory [12-14]. 

Our work is organised as follows. In the chapter 2 we introduce the action and equations 
of motion for the Einstein-Gauss-Bonnet gravity, dual to a viscous fluid in the field theory. 
Using black hole solution to the equations of motion and AdS/GFT dictionary we find an ideal 
fluid (non-dissipative) part of the energy momentum tensor. In the chapter 3 we introduce 
derivative corrections to the gravity dual corresponding to the viscous contributions in the field 
theory side. Solving the equations of motion up to second order in derivatives we reproduce 
the known result for the shear viscosity over entropy density (1.8) and show explicitly the 
breakdown of the universal relation (1.10). 


2 AdS Gauss-Bonnet black holes 

The equations of motion for Gauss-Bonnet gravity can be derived from the following action, 

5 = 5egb +‘S'b (2.1) 

where 5 egb was defined in (1.7) and Sb are appropriate counter terms, imposed to make the 
variational principle well defined [19, 20]. The Roman indices m, n = 0,... ,4 refer to bulk 
(AdS) quantities while Greek indices = 0,..., 3 refer to boundary quantities. We set 
L = 1 from now on. The equations of motion that follow from (2.1) are given by 

1 9 

Rmn - ^Rg mn - 65 mn - -j^gmuLcB 

+ 29 - 2mRmpnq - 2RjRgn + RRmn) = 0 . ( 2 . 2 ) 

The black hole solution to (2.2) is given by 


ds^ = —r‘^f{p)dt^ -|- (dx*)^ -|- 2s{p)dtdr , 


(2.3) 
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where p = br and 


s{p) = 


1 + VI - 86» 






40 


0 " 


(2.4) 


with b an integration constant associated with the mass of the black hole [ 21 , 22 ]. 

Under the gauge-gravity duality, the black hole solution is dual to a thermal state of the 
black hole. The energy momentum tensor of this state can be computed using the standard 
prescription, [14, 23, 24], 


Tnu = lini —— 
r^oo SttGs 


+ 20 — V7/i^) — + 


- 2^2 + SS 0 -I 
So 


7m. (2.5) 


where: 


7 mn = 9mn “ NmNn is the boundary metric 

Nn = S^n!\fW^ a- unit outward vector to the boundary 

= — ^ i'^niNn + n^m) IS the extrinsic curvature on the boundary (the covariant 
derivative here is with respect to the bulk metric Pmn) 

Jmn = 4 f2/C/C^„/C\ + }Cr,sJCP^}Cmn - 2/C^„/CP^/C,„ - 1C^IC„ 


K, and J are the trace (using to contract indices) of ICmn and Jmn respectively 


Inserting (2.3) into (2.5) and using the following relation 


T = 


bTT 


( 2 . 6 ) 


between the temperature of the dual theory and the horizon parameter b, derived in (A.2), 
we find 

2^4 4 

{4u>^u‘' + . (2.7) 




IGttGs^s 

The parameter G 5 is theory dependent. For instance, in the planar limit of strongly coupled 
M = 4 super Yang Mills one finds 


G5 = 


TT 


2 A ^2 ■ 


( 2 . 8 ) 


Gauss-Bonnet coupling 0 can be regarded as 0{N) corrections to the action [12, 14]. 
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3 Solving the equations of motion perturbatively 


In order to compute the transport coefficients of the fluid dual to the theory dehned by (2.1) 
we use the method of [25]; The metric (2.3) may be boosted with velocity ^1, 

resulting in the line element 


ds^ = —r‘^f{br)u^Uudx^dx^ + r'^P^ydx'^dx'' — 2s{br)u^dx^dr . 
We now promote and 6 to be a slowly-varying helds, 


(3.1) 


/3i ^ /3i(x"), b -)■ 6(x"). (3.2) 

Inserting (3.2) into (3.1) is in general, not a solution to the equations of motion (2.2). There¬ 
fore we correct (3.1) by adding new terms to the metric. The corrections to the metric may 
be classihed as follows. Let us decompose the metric into scalar, vector and tensor modes of 
the SO{3) C SO{3, 1) symmetry under which is (locally) invariant, i.e., 

ds'^ = r'^kUfj^Uydx^dx'^ + P^i^dx^ dx'^ — 2Su^dx^dr 

-|- (u^W -|- UyV^) dx^dx'^ + r'^Ii^^dx^dx’^ (3.3) 


and expand k, S, and Il^jy in gradients of and b. It is convenient to denote the n’th 
order correction to k, S, V)i and Il^i/ with a superscript (n). For instance. 


k = k^^^ -|- k^^^ + k^'^^ -|- ... = —f{p) + k^^^ -|- -|- ... 

S = + 5 ( 2 ) + ... = s{p) + + 5 ( 2 ) + ,,, 

v, = vi^^ + v^^^ + v^^^ + ... = v^^^ + vi^^ + ... 
u,, = ng) + nW + n^2) ^ ^ ^ ^ ^ + n® +..., 


where f{p) and s{p) are given in (2.4). Inserting (3.4) into (3.3) and using the equations of 
motion (2.2) we hnd, at each order in n > 0, 



Wf{p)) 


p^sl-^ef{p)) yW'j =vW 

p^fip) 


Si-^9fip) 


ttWM = p(’^) 

flU I uu "> 


fj.u 


(3.5) 
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where the quantities on the right-hand side of (3.5) depend on the known 0 to n — 1 order 
metric components. Integrating (3.5) leads to 

poo 

Si^) = - S^^\x')dx' 

Jp 

Un) ^ //kW(xOdx^ _ 2f{p) („) 

(5g - 46»/(p)) So / (5-2 _ 46»/(p)) 

vi-)=-r 

J O 


dx 


x^ (S^ - 40f(x)) Jp x^ (S^ - 46»/(x)) 


nW = - 


{So - 4(9/(x)) /f Pji) {x')dx' 


dx. 


x^f{x) 

The boundary conditions for the integral were chosen to ensnre that 

1. The boundary metric is flat: lim ds^ = r'^ppudx^dx’^, i.e. for n > 1: 

lim = 0, lim = 0, lim = 0. 

-r—-r —^ r—^ 


(3.6) 


2. The metric is regular at r = 1. 

3. is in the Landau frame: there are no 1/p^ dependence in the near-boundary ex¬ 
pansion of and with n > 1. 

With a perturbative solution to the equations of motion at hand one can use the pre¬ 
scription (2.5) to compute the energy momentum tensor order by order in the derivative 
expansion. 


3.1 First order 

Using (3.4) to first order in derivatives, we find (3.5) with 


= 0 


= b 1 

(^Sop^S^o-^Ofip))'] d. 

Yj'> = b{Sop^ (52-40/(p)))'u“d„ 

- b ( 


^ pu " ^ 

,si-4ef{p)) 


Inserting (3.7) into (3.6) we get 

5^^) = 0 

kW = 5^ 

3p 

yW = daUp 

= bTT{p)crpu , 


(3.7) 


(3.8) 
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where 


, > r So" - 4»/(i) - , 


Evaluating (2.5) by using the first order metric solution and temperature relation (A.2), 
we get 


rpfjil/ _ 


r^TT^ .. r^TT^ (1 - 80) 


IOttGsSq 

Using (1.5) we can read the shear viscosity: 




IOttGsSq 


fj 




(3.9) 


T^vr^ (1 - 80) 
IhvrGs^Q 


(3.10) 


Note that using s = ^ together with the relations (2.4), we recover the result of [13] 

2 = ^(1-SO). (3.11) 

s 47r 

3.2 Second order 

Inserting the zeroth and first order solutions (2.4) and (3.8) into (3.4) and solving (2.2), one 
finds that the solution to the equations of motion take the form 


where 


S = So + S4{p) ©4 + hip) ®5 

k = -f{br) + 6^ + ko(p) S3 + ki(p) ©i 

3p 

+ k3(/9) ©3 + k ^ ip ) ©4 + hip ) 65 

s 

Vf, = -b - daU^ + W4(p) V4^ + %(p) V5^ 

^ ^ 


+ Vi(p) +V2(p) 5J2ai + V3(p) QJ3^ 
3 

= b Trip) + ^iip) 1 

i=0 


S3 = ^P-Pdadpb, 

©1 = Vu^VUa, ©3 = idf^uh^ 
©4 = IfM^, ©5 = , 


(3.12) 


V = u^dfj,, h = e°‘^'^^Uadpu-y, 
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V4i/ = 




V 5 , 


= P>^,p^^d^dpu^, 


QJli, = ^2u = (^aP'yfj.U^'T^ul^P, ^3iy = -(TauPu° 


(i) 

and cT^jy, S^jy and a;^,y were defined in (1.4). 

The functions Si{p), ^i{p), '^i{p) and ^i{p) are determined by solving (3.6). We have not 

written the explicit form of the sources and The sources for the second order 

( 2 ) 

contribution to 11^^, P)iv ■, are given in Appendix B. 

Once the functions Sj(p), kj(p), Wj(p) and TTi{p) are determined, we can expand them 
around the asymptotic boundary, 


5 = S'o + -2 (s4 64 + S5 65) 

/c = -1 + ^ + 6^ daU°‘ + ^ (ko S3 + ki 61 + ks 63 + k4 64 + ks 65) 

p^ 3p p^ 

So b‘^ 

Vfj, = -b - U^daU^ + ^ (W4 V4^ + W5 V5^ + Vi + V2 ^3 m) 

p p^ 


P^lu = b 


TTi 7r4 

I 4 I < 7 ^ 1 / 
p p^ ' 


+'>’E 


i=0 


^ + vW 
P^ P^ ) f ’ 


where 


K 4 = 


Sq _ a TT — f9<?2 if 

, TTi — ho, vr4 — — (2 do — 1) 


252 - 1 


(3.13) 


(3.14) 


are found using the zeroth (2.4) and first (3.8) order metric functions. The coefficients vrj, 
7f2i and 7f4i can be obtained by expanding (3.6) near the boundary located at p —)• oo. In 
particular, we find that 

_ (2^g - 1) f 6Pr(t) ^ 


^4i — 


24 


t 


+ 5Pr(t) + tPp'(t) 


+ 


+ 6Pf ^"(z) + zPf ^'"(z)^ dz 


(3.15) 


( 2 ) 

where each P) Ms a part of the following decomposition 

Pg) = pg^sg) + pg^sgj + pg)sgj + pg^sg). 


( 2 ) 

The explicit form of the source terms P) (1 = 0,..., 3) appear in the Appendix B. Inserting 
(3.13) into (2.5) we find 


(3.16) 


= P (4uV + ^ AiS 

i=0 




(3.17) 








where P and rj were determined in (2.7) and (3.9) respectively and 

_ it4* {2Si - 1) 

4nb^G,So ■ 

We have not managed to reduce (3.15) to a closed form expression. We do note, however, 


that 


where 


- + APf'>-Pf = -^S{p), 


86 + Sjp^ - AOp^f jp)) 4 (gg - Sjp^ - 46/ jp)) 

Sl{l-W)p‘^{Sl-Wf{p)) So{Sl-Wf{p)) 
sep {2Si - f jp)) (5g (1 + p3) - 40/ (p)) 

Sl{l-89){Si-Aef{p)) ^ ■ 

Plugging (3.20) into (3.6) and then using (3.18), we find that 

vr^T^ ’V2 (1-^1-80 + 160)0' 


— 2Ao + 4Ai — A 2 — 


vrGs 


(1 + 


3/2 


(3.19) 


(3.20) 


as advertised in (1.10). This is the main result of our paper. 

We note in passing that the authors of [18] have found closed form expressions for the 
Ai’s. We have compared our integral expression (3.15) with their results both numerically 
and analytically and found excellent agreement. For instance, expanding (3.15) at small 0 we 
get 

2 - log 2 + (-21 + 5 log 2) 0 + ^5 + y log 2^ 0^ + 0(0^) 

1 

1-70+ —02 + 0(03) 

A 2 = [2 log 2 + (14 - 10 log 2) 0 - (28 + 21 log 2) 0^ + 0(0^)] 

= [- 1129 - 809 =+ 0 ( 0 =)]. 

with ^ ^ 

- 2Ao + 4Ai - A 2 = :^^320 0^ + 0(0^). 
ozvrGs 


“ 327rG5 
TT^r^ 
“ 327rG5 
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Appendix A Hawking black hole temperatnre 


To find the black-hole temperature we convert our ansatz (2.3) into a non-zero temperature 
black 3-brane metric. With dt ^ dt — one gets 

dr^ ■ 9 

ds‘^ = -r'^f{br)dt^ + „ -k [dx'-) . (A.l) 

r^f(br)/S^ 

Since the Hawking temperature is completely determined at the vicinity of horizon, we expand 
the metric functions near r = (in our case rh = with /{br^) = 0) 


c2 J 2 

= - ir^fibr)y\ (r - rh)dt‘^ + ^ ° -r 

tr-rn (r^/(6r))'| (r - r/,) 


+ r^ {dxy 


Defining r — 


(r^f{br))'\ 


■4^ 


^ and using the zeroth order solution (2.4), we find 

4aV252 


ds = — 


° z^dt^ + dz^ + {dx^y . 


62 


After a Wick rotation to Euclidean time r, we have (without the last term) a flat space metric 
in cylindrical coordinates ds^ = dz"^ + z'^dip'^ where (p = To avoid a conical singularity 

at z = 0, ip must have a periodicity of 27r. Since periodicity of Euclidean time is the inverse 
temperature, we have 27r = ^ and therefore 


T = 


So 

bTT 


(A.2) 


Appendix B Second order sonrce terms 

(‘2'\ 

Source terms for the second order IVpp metric functions, see (3.16) and (3.6). 


>{ 2 ) 


S^ (24 0 + (1 - 80)/) S^ (-40 0 - 3(1 - 80)p^) 


/ (S^-4df(p)y 


+ 


p2 (ss-40f(p)y 


7r(p) + 


-2Sop^ 


s^-mp) 


Tr'(p) 


p(2) _ P 

^ S^\l-80 


-40 ^ -25g(l-80) ^ 


354 


^ -55(400 + 3(1 - 80)/) ^ ^ 
H-—W-- -T^{p) 


(52 - 40/(/)' - 40/(/; p2 (52 - 40/(/) 


+ 


+ 


2 / 


-2 -5S (1 - 80) 

+ , „ -/ + 


352 


tt' 


TT (p) 


+ 


So \^l-80 (52-40/(/)' 5o2-40/(/y 

/ (-165^0 + f{p) (5g - 40/(/ (5 o 2 (1 - 160) + 40/(/) (35o2 - 80/(/))) 

52 {1-89) (Sl-4()f{p)f 

. , -49p«/(p) (2S„= - f(p)) „ , , 

^ (p) + +271-7777+2-++++7’^ (0)^ (0) 


y'ip))‘ 


-1 


+ 


^0^ 


5o Vl-80 5 o2-40/(/ 


52(1 -80) (52 -40/(/) 
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2SIp , 2S^ (40 0 + 3(1 - 80)p4) ^ , 4 ^ 00 ' , 

s§-4nm ^ ^ ss-mp}" 


(2) ^ 5g(240 + (l-80y) SI (1 - 80) (35g - 160)p ASj (69 - Sl)p 

3 0(i-80)p3 e{sl-4ef{p)f o{sl-4ef{p)) 
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